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Abstract. A deposited drop of bovine serum albumin salt solution experiences both gelation and fractura- 
tion during evaporation. The cracks appearing at the edge of the gelling drop are regularly spaced, due to 
the competition between the evaporation-induced and relaxation-induced stress evolution. Subsequently, 
the mean crack spacing evolves in an unexpected way, being inversely proportional instead of proportional 
to the deposit thickness. This evolution has been ascribed to the change with time of the average shrinkage 
stress, the crack patterning being purely elastic instead of evaporation-controlled. 

PACS. 6 2.20.Mk, 47.54.+r, 82.70.Gg 



1 Introduction 

The study of crack dynamics in drying materials meets 
the basic interest in nonequilibrium pattern formation and 
the industrial need of controlling drying processes (paint 
technology, . . . ). But whereas the understanding of crack 
behaviours is growing, partly due from the interplay be- 
tween mechanics and statistical physics approaches pQ, 
the investigation of the cracking of drying complex fluids 
is just setting in. 

As part of a general study of biocompatibility between 
prosthetic materials and real biological liquids, we focussed 
on the evaporation of a protein solution drop and on the 
evolution of the solid deposit. Although the gelation pro- 
cess of our biological liquid and the initial stage of the 
. crack pattern formation show similarities with like phe- 
nomena in mineral colloids |2lri| , the characteristics of our 
experiments yield further insight in the following steps of 
the cracking. 

Our working liquid is a phosphate buffered saline so- 
lution (pH = 7.4, ionic strength / = 0.2 M) of bovine 
serum albumin (BSA) at various concentrations, which 
may be classified as a colloidal sol with the BSA molecules 
amounting to ellipsoidal particles of dimensions 4 x 4 x 14 
nm 3 4 . The volume fraction of protein is always a few 
percents in our concentration range. In the course of the 
experiment, a 15 mm 3 sessile drop is deposited onto a glass 
slide. Great care is taken for the cleaning of the substrate, 
residual deposits being liable to disturb the wetting and 
cracking behaviour. Keeping in the surface tension-driven 
regime, tensile drops have been experienced to exhibit the 
same behaviour. Immediately, proteins are adsorbed on 
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the substrate, leading to a strong anchoring of the triple 
line 0. 

The first step of dessication is then characterized by 
a fixed base radius of the drop, a constant evaporation 
rate all over its free surface and a spherical cap shape in- 
duced by surface tension (Bond number smaller than one). 
These three features have been shown to remain compat- 
ible thanks to outward flows inside the drop 0- Subse- 
quently an apparently solid ring-like "foot" progressively 
forms close to the edge of the drop, and the three phase 
line migrates inwards, provoking the extension of the solid 
until complete dessication. 

Due to the competition between the adhesion of the 
deposit on the substrate and its evaporation-induced or- 
thoradial shrinkage, the foot rapidly experiences regularly- 
spaced radial cracks (see Figure 0. The periodicity of 
the fractures in such evaporating colloidal solutions has 
turned out to originate in the interplay between the stress 
relaxation induced by the crack and the stress enhance- 
ment yielded by the water drainage from the inner liquid 
to the crack surface . 

But unlike the case of mineral colloids, the subsequent 
evolution of our crack pattern is characterized by the stop 
of a growing number of fractures, always preceded by a 
right angle turn. In the same time, purely orthoradial 
cracks appear (see Figure |2J. Finally all radial fractures 
cease to propagate simultaneously and, depending on the 
protein concentration, two scenarii may occur in the cen- 
tral zone: Small-scaled disordered fracturation or dendritic 
crystallisation (see FigureOJand 8J- Respective behaviours 
as a function of the BSA concentration are gathered in Ta- 
ble □ 
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2 Dynamics 

Before studying the fracturation, we concentrate on the 
evaporation dynamics. To evaluate our dessication time 
to = ^ with ve — "■jTjf water flux per unit sur- 
face (i?o initial base radius, S surface and V volume of 
the drop), we have performed measurements of the evo- 
lution of the drying drop height h with time t (see Fig- 
ure . Considering a constant base radius and a spher- 
ical cap shape and choosing for S the overall free sur- 
face of the drop, including the foot where evaporation still 
takes place 0, the dessication time can simply be written 



Furthermore, weighing measurements shown in Figure 
have brought an additional evaluation of the dessica- 
tion time with the above-mentionned assumptions: rrj = 

— ,d™"f° with the initial air-drop surface, m the mass 
y — )t=° 

of the drop and p the density of the solution, determined 
by picnometry. The obtained numerical value tjj = 1.5 x 
10 4 s fits with the previous determination. As can be 
stated in Figure^! evaporation dynamics compares tightly 
in pure water and protein solution. Indeed, the volume 
fraction of the protein keeping small, td is quite indepen- 
dent of the BSA concentration. 

Then we have focussed on the properties of the solid 
deposit. Such protein solutions are known to gel during 
dessication, helped by the increasing ionic strength To 
check the physical nature of the obtained biopolymer gel, 
we placed the slide in water and observed complete redi- 
lution of the deposit. Measurements of adsorption spectra 
and surface tension of the sol before and after the gelation 
and redilution process give identical results, proof of the 
non-chemical nature of the bonds in the solid. 

To evidence the sol- gel transition kinetics, we have car- 
ried out electrical conductivity measurements of the evap- 
orating drop by placing electrodes on the slide, at two 
edge points of the drop. As can be seen on Figure the 
progressive hindering of ionic conduction through the liq- 
uid by the connection of the protein molecules gives rise 
to a decrease in the current. We have taken as gelation 
time tq the elapsed time at the achievement of the low- 
conductivity plateau, tq is found to decrease smoothly 
when concentration rises (see Table ^1 . 

Finally the cracking characteristic time at which the 
first fracture is detected has been measured and found 
quite independent of concentration: tq — 8 x 10 2 s. 



3 Initial cracking 

We now point out the fracturing behaviour of the gel. The 
fact that the central part of the drop remains liquid and 
the three phase line anchored while the foot starts shrink- 
ing results in orthoradial stress, then in radial cracks. In 
a quantitative theory of the ID fracturing behaviour of a 
mineral colloidal suspension confined in a Hele-Shaw cell, 
Limat & Allain have found that the crack spacing A should 
be roughly proportional to the cell thickness |7| . We have 



implemented this model for our open geometry, consider- 
ing that the same stress field exists in our drop foot as in 
one longitudinal half of their slab, leading to an equiva- 
lent cell thickness of twice the foot height. The theoretical 
and experimental initial crack spacings Ao are reported 
in Table ^ Considering the degree of approximation of 
the model and the difficulty for the physical parameters 
evaluation, agreement is found satisfactory. 

The foot heights listed in Table have been measured 
with a holographic interferometry device, thanks to an 
original geometrical analysis ^UJ- With this technique, a 
reference hologram of the studied surface is first recorded 
with the use of two beams (reference and object), the in- 
cident angle of the object beam being a. This hologram 
contains information on the amplitude and phase of the 
reflected object beam, enabling a 3D reconstruction of the 
surface. Then a mirror is tilted, changing the incident an- 
gle of the object beam to a + Aa and a new hologram is 
recorded. The interference of the two holograms results in 
a fringe pattern, an example of which is shown in Figurc0 
We have shown in ^J] that flat surfaces lead to fringes par- 
allel to the tilt axis, whereas any defect of height m in the 
surface creates a fringe shift p with p = mx tan (a + Aa). 
The deposit topography has therefore been derived from 
the measurement of the fringes deviation on the interfer- 
ence pattern caused by a tilt Aa = 0.5° from the incident 
angle a = 80°. Validation of the results have been pro- 
vided by complementary profilometer measurements. This 
method has moreover permitted to check that the flatness 
of the slide was better than 2 /urn over the distances in- 
vestigated here. 

4 Cracking evolution 

After the initial breaking, number of cracks merge during 
propagation and the mean crack spacing A may increase. 
In the vicinity of a crack, the stress is purely parallel to 
its edges, so when a radial fracture tends to stop, it expe- 
riences a turn to join an adjacent one with a right angle. 
The evolution of the mean crack spacing A (taken as the 
ratio of a perimeter to the number of fractures crossing it) 
when going from the center to the edge of the dried drop 
is shown in Figures |H1 and El with same scales for 40 and 
60 g.l -1 BSA concentrations. 

As has been qualitatively JJ| and quantitatively (see 
above) demonstrated in different configurations, the crack 
spacing should be roughly proportional to the layer thick- 
ness. We have therefore used our holointerferometry mea- 
surements to get informations about the thickness of the 
deposit along the radius of the dried drop. Results are rep- 
resented in Figures|H|for the 40 g.l -1 case andEDfor the 60 
g.l -1 case and give a clear view of the fact that the link 
between crack vanishing and layer thickness is different 
from what could be expected. 

In fact, if we evaluate the characteristic time of dif- 
fusion of water in the gel between two cracks during the 
initial stage, the order of magnitude is found greater than 
10 minutes, i.e., comparable to the whole fracturation 
time (~15 min.). So after the initial cracking, evaporation 
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ceases to constitute a leading feature of the phenomenon. 
Therefore we must now look back at the involved kinetics 
to assess the governing mechanisms. 



5 Crack pattern 

For an identical time hierarchy (tq < tq < td), Pauchard 
et al. observed disordered crack patterns during the evap- 
oration of an aqueous silica sol deposited drop, caused by a 
buckling instability 0. This one consists of a distortion of 
the drop shape caused by the gelation of the surface when 
the inner sol is still liquid. We do not face this morpholog- 
ical behaviour, the central zone solidifying progressively 
through the propagation of the gel front. 

Before ascribing this discrepancy to the respective na- 
tures of the mineral and biopolymer gels, the influence of 
the geometrical parameters has to be investigated. Show- 
ing an axial symmetry, only two dimensions have to be 
probed. First the drop volume has been varied from 2 to 
25 mm 3 with solutions of concentrations 20 to 60 g.l -1 . 
Then, for identical concentrations, the contact angle has 
been changed by testing several substrates of biomedical 
interest, the wettability characteristics of which had al- 
ready been studied 5 : nitrogen implanted (59°) and vir- 
gin (62°) titanium based alloy TA6V, GUR 415 (65°) and 
Chirulen (67°) ultra high molecular weight polyethylenes, 
sintered alumina (66°) and zircona (67°), stainless steel 
316L (68°), and a chromium-cobalt alloy (71°). Addition- 
ally, silane-coated glass (>90°) has been tested. Whereas 
the fractures penetration depth variation with concentra- 
tion and contact angle show interesting features — which 
will be studied elsewhere — , in all cases the drop morphol- 
ogy and cracking evolution are similar to the case of the 
15 mm 3 drop deposited on glass. One only difference can 
be noticed for volumes smaller or equal to 5 mm 3 and 
silane-coated glass where all initial radial fractures stop 
propagating at the same time instead of hierarchically due 
to the small base radius of the drop. 

So the universal radial nature of the crack pattern (in- 
stead of a disordered nature for a mineral gel) seems to be 
independent of geometrical considerations and can reason- 
ably be attributed to the physical nature of the biopoly- 
mer gel, the latter exhibiting smoother mechanical prop- 
erties than chemical gels. Its low tensile strength in the 
gelling central zone should thereby result in an overall 
stress in the gel front keeping an orthoradial direction. 



6 Mean crack spacing evolution 

The dessication kinetics being non-critical after the ini- 
tial cracking, information is now transmitted elastically 
between cracks, through the flux of displacement field, as 
pointed out in impact fracture experiments |12j . which is 
a much faster mode. Thus we face a pattern selection sys- 
tem comparable to the one encountered in thermal crack 
growth experiments J3| : The release of elastic energy dur- 
ing the propagation controls the crack spacing. 



Indeed for the propagation of a crack to be possible, 
the elastic energy E e liberated by the breaking must at 
least equal the surface energy E s required by the two 
formed lips (Griffith's criterion). We are considering here 
as a first approximation purely brittle breaking, dissipa- 
tive effects (plasticity, friction) being neglected. During 
the progression 51 of a fracture, E e is proportional to the 
relaxed volume e x A x 51 and E s to the created surface 
area e x 51. So the fulfilment of the E e = E s condition 
does not imply, in the case of a change of e, a modifi- 
cation of the crack spacing. Therefore the evolution of A 
when approaching the center of the drop, if any, should be 
related to a change in the only other varying factor along 
the fracture path: the stress field. 

For example, in the 40 g.l -1 case, A evolves roughly 
as r~2 (r distance from the center of the drop), so linear 
elasticity yields a change of the stress a proportional to 
A~ 2 so to r 4 . This slow decrease of a when going from the 
edge to the center of the drop has to be related to the de- 
crease of the average stress with time: Cracks meet during 
their propagation a gel more and more relaxed by the con- 
tinuous fracturation |14j . In the 60 g.l -1 case, A exhibits 
no actual evolution, which should stem from a constant 
stress during the cracking process. This constancy of the 
stress is likely to be a consequence of the large and quite 
uniform layer thickness, which hinders relief mechanisms. 

So the dynamics of the mean crack spacing has re- 
vealed itself to be actually related to the deposit thick- 
ness, but not in the expected way: It is rather inversely 
proportional instead of proportional. Indeed, dessication 
playing no more role, only cracking influences the crack- 
ing geometry through stress relief, which is made easier 
by decreasing layer thicknesses. 



7 Conclusion 

We have investigated the evaporating, gelling and crack- 
ing behaviour of a deposited drop of protein solution. 
The initial stage is characterized by the appearance of 
regularly-spaced cracks at the edge of the gelling drop. 
This ordered patterning was known to originate in the 
interplay between the water diffusion towards the crack 
surface and the elastic relaxation due to the fracture. The 
general radial patterning was found to be a corollary of 
the soft nature of the BSA gel. Then a transition from 
this evaporation-controlled behaviour to a purely elastic 
one occurs, where the change of the crack spacing has 
turned out to be a consequence of the shrinkage stress 
evolution with time. 

This first study of the fracturing behaviour of a drying 
biological liquid has shown similarities with other com- 
plex fluids. But discrepancies have also been emphasized 
and an investigation of the specific mechanical properties 
of the biopolymer gel is still needed to confirm the as- 
sumptions made on the dynamics of the fracturation. The 
influence of the substrate on the cracking, especially on 
fracture penetration depths, constitutes also a promising 
way of investigation. 
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Fig. 1. Beginning of the foot fracturation of a drying 40 g.l 1 Fig. 2. Following of the fracturation of the drying drop of 
solution drop 840 s after deposit. Figure [I] 500 s later. 
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Fig. 3. Dendritic crystallisation of the central zone of a dried 
10 g.l -1 solution drop. 
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Fig. 4. Time evolution of a drying 40 g.l - solution drop 
height. As a matter of comparison the values for a pure water 
drop are also displayed. 
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Fig. 5. Time evolution of a drying 40 g.l -1 solution drop mass. 
Values for pure water are also displayed. 
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Fig. 8. Layer thickness e and crack spacing A as a function of 
the reduced distance r/Ro from the center of a dried 40 g.l -1 
drop. 
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Fig. 6. Time evolution of the electric current between two edge 
points of the drying 40 g.l~* solution drop. 
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Fig. 9. Idem as Figure |H| (same scales) for a dried 60 
drop. 
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Fig. 7. Holointerferogram of a dried 40 g.l 1 drop deposit. 



(g-r 1 ) 


foot height 

(H 


TG 

(10 3 s) 


■» exp 

A o 


\ theo 
W 


final behaviour 


20 


28 


3.8 


100 


124 


crystallisation 


30 


31 


3.7 


120 


135 


crystallisation 


40 


36 


3.6 


160 


154 


fracturation 


60 


81 


3.4 


250 


309 


fracturation 



Table 1. Foot height, gelation characteristic time, experimen- 
tal value of the initial crack spacing, theoretical value of the 
initial crack spacing, final behaviour of the deposit center as a 
function of BSA concentration. 



